Copyright © 1976 American Telephone and Telegraph Company 

The Bell System Technical Journal 

Vol. 55, No. 4, April 1976 

Printed in U.S.A. 



Analysis of a Gradient Algorithm for 

Simultaneous Passband Equalization 

and Carrier Phase Recovery 

By D. D. FALCONER 

(Manuscript received December 11, 1975) 

A two-dimensional receiver structure has been proposed, incorporating 
two innovations: passband equalization, which mitigates intersymbol 
interference, and data-directed carrier recovery and demodulation following 
equalization, which enables compensation of carrier frequency offset and 
phase jitter, but does not require transmission of a separate pilot tone with 
the data signal. The receiver is fully adaptive; the adjustment of the equal- 
izer tap coefficients and of the estimate of the current channel phase shift is 
based on a gradient algorithm for jointly minimizing the mean squared 
error with respect to those parameters. 

In this paper, we analyze the dynamic behavior of the deterministic 
gradient algorithm {where channel parameters entering into the gradient 
expression are assumed known in advance). The corresponding estimated 
gradient algorithm (where these parameters are initially unknown) has 
previously been studied experimentally, but is not treated here. 

The first part of the present study concerns system start-up (or transient) 
response when the channel's phase shift is fixed. Examination of the analyt- 
ical solution leads to the qualitative conclusion that, if the equalizer tap 
adaptation coefficient /3 is small relative to the phase-tracking coefficient a, 
the added phase estimation feature does not strongly affect the start-up 
behavior of the passband equalizer under typical operating conditions. 
Indeed, if the equalizer tap coefficients all start at zero, their evolution in 
the deterministic gradient algorithm is completely unaffected by the phase- 
tracking loop. 

The second situation analyzed is the steady-state response of the system 
to a constant carrier frequency offset. In this case, the phase-tracking loop 
is found to reduce the resulting rate of rotation of the equalizer taps to about 
0/(a + 0) of the original frequency offset. As a result, the degradation in 
system mean squared error due to frequency offset is typically quite small. 

The final analysis is of the response of a linearized version of the 
receiver structure to sinusoidal phase jitter. When the channel's linear 
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distortion is not too severe and the coefficient is small, the system mean 
squared error owing to phase tracking error is found to approximate that 
of a simple, first-order, phase-locked loop. 

I. INTRODUCTION 

The combination of adaptive equalization and decision-directed 
estimation of a fixed carrier phase offset in suppressed-carrier pam 
modems by means of a gradient algorithm has been suggested by 
Chang 1 and by Kobayashi, 2 the latter also including adaptive timing 
recovery. The receivers contemplated in those papers demodulated the 
received data signal prior to equalization and carrier phase estimation. 

Reference 3 describes an alternative receiver configuration for two- 
dimensional modulated data transmission systems, combining equaliza- 
tion and carrier recovery. This receiver's distinction is that it employs 
a passband equalizer 4 whose reference signal consists of receiver 
decisions amplitude-modulating a carrier whose phase shift is the 
receiver's estimate of the channel phase shift. Following the passband 
equalizer is a demodulator which compensates for the channel's phase 
shift (which may be time-varying as a result of frequency offset or 
phase jitter). 

The receiver's estimation of the carrier phase shift is based on a 
decision-directed gradient algorithm for estimating a fixed phase shift, 
as proposed in Refs. 1, 2, 5, and 6. An advantage of the demodulator 
following the equalizer is that the demodulator's phase reference is 
delayed relative to the actual channel phase shift by only one symbol 
interval instead of by the entire equalizer delay as in the traditional 
"baseband" receiver configuration. This fact, plus the provision of a 
sufficiently large gain coefficient in the phase-tracking gradient algo- 
rithm, makes possible tracking and compensation of typical conditions 
of frequency offset and phase jitter that may occur on voiceband 
telephone channels. Computer simulations, reported in Refs. 3 and 7, 
have confirmed this capability. 

In this paper, we study the dynamic behavior of the gradient algo- 
rithm for jointly adjusting the equalizer tap coefficients and the phase 
estimate in each of the following situations: (i) start-up (transient 
response) for a fixed carrier phase shift ; (ii) steady-state response to a 
frequency offset; (Hi) steady-state response to sinusoidal phase jitter. 
Throughout, we consider only the deterministic gradient algorithm; 
that is, receiver decisions are assumed perfect, and the gradient of the 
mean squared error as a function of equalizer tap coefficients and carrier 
reference is assumed known. A stochastic gradient algorithm, which 
would be used in practice, has been simulated, 3,7 but is not treated in 
this paper. 
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II. SYSTEM EQUATIONS 

The transmitted two-dimensional modulated data signal is assumed 
to be of the form 

«(*) = Re {£ A n g(t - n7> 2 "«'}, 



where A„ is a two-dimensional (complex-valued) data symbol trans- 
mitted in the nth symbol interval, g(t) is a band-limited baseband 
pulse waveform, T is the duration of a symbol interval, and f c is the 
carrier frequency. The set of possible discrete complex values that each 
An can assume constitutes the signal constellation. Quadrature ampli- 
tude modulation (qam) and digital phase modulation (pm) systems are 
familiar examples of two-dimensional modulation systems. We shall 
assume that successive data symbols are uncorrelated ; i.e., 

(A„A* m ) =1 f or n = m 
= otherwise. 

Figure 1 shows the receiver structure. The received signal, after 
transmission through a noisy, dispersive channel which may introduce 
a slowly time-varying phase shift, is passed through a phase splitter 
to produce parallel in-phase and quadrature components. These parallel 
waveforms can be represented as a single complex waveform that is 
sampled and passed on to a passband transversal equalizer with, say, 
2N + 1 complex-valued tap coefficients. In the nth symbol interval, 
when a decision is to be made on the nth data symbol, the latest 
(2N + 1) complex- valued samples stored in the (2N + l)-tap pass- 
band equalizer can be represented by the complex (2N + l)-dimen- 
sional vector R n e jSn , where 9„ is the channel phase shift (assumed 
quasi-stationary in the nth symbol interval). A sequence {0„} changing 
at a constant rate with time is an example of frequency offset, while 
{On} varying randomly or quasi-periodically constitutes phase jitter. 
Typically, the change in 6 n in one or two symbol intervals is so small 
as to allow us to neglect the phase-to-amplitude modulation con- 
version effected by filtering the sequence of incidental frequency- 
modulated components [e' Bn \. 
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Fig. 1 — Two-dimensional receiver. 
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The (2iV + 1) complex equalizer tap coefficients in the nth symbol 
interval are denoted by the complex (2N + 1) -dimensional vector 
C„ = (c-n, ■ ■ -, Co, • • •, cat).* The symbol * will denote transposed 
complex conjugate throughout. Then the nth complex equalizer output 

Q n = C:R n e*«, (1) 

the real part being interpreted as the in-phase component and the 
imaginary part as the quadrature component. 

The receiver's estimate of 6 „ is a real quantity denoted by 0„, and 
the demodulator output is written 

Y n = Q n e-M*'° n ™»\ (2) 

This quantity is passed into a simple quantizer to produce A n , which 
is the receiver's decision on A„. Based on this decision, the complex 
reference signal used for updating the equalizer taps and the phase 

estimate is 

Q n = i n e* 2 "<»™»>. (3) 

We define the properties of the channel in terms of expectations 
(denoted by < )) with respect to the ensembles of information symbol 
sequences and additive noise samples. The complex impulse response 
X is defined by 

X = jr^-r* (A:R n )e-*"<»T. (4) 

The positive definite Hermitian a matrix of the channel is defined by 

a = (jx^> < R - R »>- (5) 

The normalized mean squared error in the nth symbol interval is 
defined to be 

<- = /TT-m < I One-* 2 ""*™ - A „ | 2 ), (6a) 

which, by virtue of (1), (4), and (5), can be rewritten as 

e n = 1 - x*a- i x + 7 n, (6b) 

where y n = ElaE n ^ is the excess mean squared error and E„ is a 
tap-error vector, 

E„ ■ C«8*k-« - CT l X. (7) 

Since a is positive definite, the value of e„ is a positive minimum, 1, 
1 — X*« -1 X, when the equalizer taps C„ and phase shift estimate 0„ 

f The positive quantity X*a -1 X is therefore less than unity/a fact which is exploited 
in the appendix. 
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are adjusted so that E„ = 0, or 

C„e'*- = Ct-^Xe*-. (8) 

This equation is also the condition for the gradients of e„ with respect 
to C„ and 6„ to be jointly zero; it is satisfied by an infinitude of points 

(C n , 0„). 

Thus, a gradient algorithm can be used to adjust the tap coefficients 
C» and phase estimate 0„ recursively toward optimal values. The 
equations governing the evolution of {C„} and {0„} are 3 

C n+1 = (/ - /3Ct)C„ + (Or**** (9) 

and 

t n+l = fi„ + alm [C:X e -^-*«>], (10) 

where / is the identity matrix and and a are positive gain coefficients. 
These equations [or the equivalent equations (13) and (14)] form the 
basis for the results in this paper. 

In practice, X and a would generally not be known in advance, and 
the following stochastic gradient algorithm, 3 involving^ the equalizer in- 
puts R„e*», outputs Q„, and modulated decisions Q n , would replace 
the deterministic gradient algorithm described by eqs. (9) and (10). 

C n+1 = C n - /3R„e#»(Q« - &)■ (11) 

n+ i = 0„ + T -£ Ta Im(Q B &). (12) 

These are coupled stochastic difference equations, since successive 
vectors }R n } are correlated random variables. Simple stochastic 
gradient algorithms have been studied by Widrow. 8 The application 
to equalizer adaptation, where no phase recovery is involved and under 
the assumption that the {R n } are uncorrelated, has been studied by 
Ungerboeck, 9 by Gersho, 10 and by Gitlin, Mazo, and Taylor. 11 The 
extension to correlated vectors (R„) has been introduced by Daniell. 12 

That the algorithm specified by (11) and (12) converges and can 
perform satisfactorily is confirmed by the computer simulations re- 
ported in Refs. 3 and 7. Analysis of the stochastic gradient algorithm is 
complicated by the possibility of a cycle-slipping phenomenon as in 
phase-lock loop systems. References 5 and 6 deal with continuous- 
time, decision-directed, phase-locked loops in the absence of adaptive 
equalization. 

However, insight can be gained by studying instead the deterministic 
gradient algorithm of (9) and (10), since the estimated gradient algo- 
rithm can be interpreted as implicitly performing the averaging in- 
volved in determining X and a, provided the signal-to-noise ratio is 
high and the gain coefficients a and /3 are sufficiently small. 
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Using definition (7), we can rewrite the coupled difference equations 
as 

E n+ i = (I - 0a)Ene'' (S » +1 - A " +l) + a-^e^n+i-An+D _ i) (13) 

and 

A n+1 = «Im(E:X), (14) 

where 

A n +i = 6 n +i — On and A„ + i = d n +i — dn. 

III. SYSTEM START-UP WITH FIXED CHANNEL PHASE SHIFT 

In this section, we study the behavior of the deterministic gradient 
algorithm during start-up, assuming the channel's phase shift is fixed : 
d n = 0. f General theorems tell us that, if the initial error and the 
coefficient of the gradient algorithm are small enough, convergence is 
guaranteed. 13 However, we are interested in sharper results for the 
specific problem at hand. 

The solution of (13) and (14) will depend on the initial choice of 
E (or Co) and O . It is interesting to consider first the special case 
Co = 0, the all-zero vector; i.e., E = — Qr l X. In this case, 

Ai = -a Im [X*a- l X] = 0, 

since Q, is Hermitian, and 

E x = -(/ - ptya-'X. 

Continuing, it is easy to show that 

A„ = for all n 

and that 

E n = -(I - payor 1 *. (15) 

Thus, at least for this special all-zero starting condition, the estimated 
carrier phase shift 0„ does not change at all and the start-up behavior 
of the deterministic algorithm is exactly the same as that of the pass- 
band equalizer alone. 4 

Let us now consider the more general case, when Eo is not necessarily 
equal to the right-hand side of (15) for some n ^ 0. We remark that 
the mathematical formulation of this start-up situation will be basi- 
cally the same as that of a system transient caused by an abrupt change 
in the channel's carrier phase shift. 

Expression (6b) for the normalized mean squared error involves the 
positive definite quadratic form EjJ(iE„ = y„- We can bound this term 



* There is no loss of generality in assuming a fixed phase shift of zero, since any 
nonzero fixed phase-shift factor e' 1 can be incorporated in the complex channel im- 
pulse response X. 
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and study its evolution by writing down a recursive expression for it 
and upper-bounding the right-hand side of that expression. Using (13) 
with A„+i = for all n, we can write 

7n+1 = E*(/ - pa)a(i - /sa)E„ + x*a- l x\e^^ - i| 2 

+ 2 Re {E:(7 - /3Ct)X(l - (r*»»)). (16) 

The right-hand side of expression (16) is upper-bounded in the ap- 
pendix. The derivation of the bound requires the following assumptions 
about the channel and algorithm parameters. 

Assumption (1) : The initial value 70 = EoCtEo is less than unity. This 
condition is fulfilled, for example, if Co = 0; i.e., E = — G _1 X, for 
then 70 = X*(i _1 X ^ 1, since the positive quadratic form X*d -1 X, 
which is one minus the minimum mean squared error, must be less 
than unity. 
Assumption (2) : a < cto, where ao is the solution of 

«o(l + V70) = 2 sine (aoVTo), 

where 

sin 6 
sine 6 = — r— • 
a 

Assumption (S) : Let the maximum and minimum eigenvalues of the 
positive definite Hermitian matrix a* be denoted respectively by X max 
and X min . Then the gain coefficient /3 must satisfy 

< fi < 2Xmi " 



X m ax(l + «o) 

where eo is defined in terms of a by 

a 1 1 + V70 + -2 ) = 2 sine (aV-yo), a < a . 

Figure 2 illustrates the solution of the equations defining el and ao- 
For example, if we assume a = 0.5 and 70 = 1, then a is 0.88 and 
el is 0.543. 

The upper bound obtained in the appendix is 

Tn+ i = E: +1 aE n+1 ^ E:ctE„ - 2/3E:a 2 E n 

+ /3 2 (1 + el)E* n a 3 E n . (17a) 

An explicit bound on 7„+i is obtained by first weakening (17a) using 
(41) of the appendix to obtain 

7«+i ^ (1 - 2/3X min + /3 2 (1 + ag)XLx)7n, (17b) 

so that 

7n+i ^ (1 ~ 2/3X min -I- 2 (1 + eoOXLxWo- (l?c) 
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Fig. 2 — Illustration of the definitions of e£ and a . 



In the absence of phase tracking, a = A„ = 0, and the mean 
squared error at step n + 1 of the deterministic gradient algorithm is 
obtained directly from expression (15) 9,10,14 as 



7n+l = EMI - p\i) 2n \8oi\ 



(18) 



where the summation is over all the eigenvalues of the matrix ft, the 
{X,} are the set of eigenvalues, and S oi is the inner product of Eo with 
the normalized ith. eigenvector. 

Comparison of the upper bound (17c) for the joint equalizing and 
phase-tracking receiver and the exact expression (18) for the equalizer 
alone yields some insight into the penalty in convergence rate imposed 
by the additional phase-tracking algorithm. Consider an example where 
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all {A,} (and therefore X max and X min ) are equal to a common value X. 
This would represent the case of a channel with delay distortion but 
not amplitude distortion (flat Nyquist equivalent frequency charac- 
teristic). Then inequality (17c) becomes 

7n + i ^ [(1 - W + /3 2 X 2 eg>yo, (19) 

and, recognizing that 

7o = JL X,|£ .| 2 , 
t 

we can write equality (18) for the case of no-phase tracking as 

7m-i = (1 - W*7o. (20) 

In practice, the equalizer adaptation coefficient is small (/3 « 1/X), 
to minimize the mean squared error resulting from a practical stochastic 
gradient algorithm. 9 Thus the right-hand sides of (19) and (20) should 
be nearly equal, and we conclude that an ideal gradient algorithm for 
joint phase tracking and equalization should not converge appreciably 
slower than the equalizer adjustment algorithm alone. An exact 
analytical evaluation of the effect of phase tracking on the convergence 
of a practical stochastic gradient algorithm for a severely distorted 
(X mft x » X min ) channel remains elusive. However, the results of this 
section suggest that the influence of the phase-tracking parameter a 
in the convergence is relatively small. This conjecture is bolstered by 
the experimental results summarized in Figs. 3a and 3b. A 9600-b/s 
two-dimensional data transmission system was simulated, employing 
the stochastic gradient algorithm described by eqs. (11) and (12). 
The transmission channel, whose frequency characteristics are shown 
in Fig. 3a, was regarded as severely distorted (it violates the minimum 
standard for private line voiceband channel data transmission). The 
plots of measured mean squared error versus time for a = and for 
a = 0.2 shown in Fig. 3b are very similar, indicating that little penalty 
in convergence rate is to be ascribed to the use of joint decision- 
directed phase tracking. 

IV. CASE OF FREQUENCY OFFSET 

In this section, we study the behavior of the system in the presence 
of frequency offset by obtaining steady-state solutions to eqs. (13) and 
(14) when the channel phase shift increases linearly with time; i.e., 
A„ = 2irAT, where A is the frequency offset. In this case, eq. (13) 
becomes 

En+1 = (/ - /3a)E n e>< 3 » + '- 2 ' AT) + a- 1 X(e^- + r- 8 * Ar) - l). (21) 
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Fig. 3a — Frequency characteristics of the simulated channel. 

A steady-state solution to (21) and (14) is obtained by substituting 
the trial solution, 

E„ = E 

An = 2tt(A + 8)T, 

and then solving for the fixed quantities E and 8. The substitution 
results in 

E = ( e > WT - l)M- l Or% (22) 

where M is the matrix 



and 



M = I - e>* rfr (I - 0a) 

2t(A + b)T = aim (E*X). 

= a Im l(r*«* - l)X*a- l M*-*X]. (23) 



It is clear from the definition of M that the eigenvectors {y,-}'^ of &, 
which form a complete orthonormal set, are also those of M. Thus, 
expressing the vector X as a linear combination of y,-, we write 

X= £ G iVi) 

i=-N 
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Fig. 3b — Convergence with and without phase tracking (ideal reference; all 
equalizer top coefficients start at zero). 



we can rewrite (23) after a little algebra as 



iV 



2tt(A + B)T = a lml(er^ T - 1)] £ 



I Gil 



= -a/3sin27r6T £ 



iJ± N X,[l - e-»* ST (l - /3X,)] ' 



-jv 1 - 2(1 - /3X.) cos 2t8T + (1 - /3X,) 2 ' 



(24) 



where {\{}«L-jv are the eigenvalues of Q. and are positive and real. 
The excess mean squared error is similarly given by 



7n = E:aE n = \e*« T - i| 2 x*a- 1 M*aM- 1 a- 1 x 

= 2(1 - cos27rfir) 

|G,-| 2 



-v 
t ]k N X,[l - 2(1 - /3X.) cos 2t8T + (1 - 0X,) 2 ] 



(25) 



Equation (24) is a transcendental equation whose solution 5 is 
clearly not zero in general. The quantity 5 may be interpreted as a bias 
in the receiver's estimate of the frequency offset. This "residual" 
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frequency offset then must be compensated for by a rotation of the 
equalizer complex tap coefficients at rate 8 Hz. 

For purposes of illustration, we again consider only a special case 
of a "good" channel, for which all X, = 1 and £,• | (?,| 2 = 1. Then (24) 
becomes 

._ —a/3 sin 27r5!T , oex 

2 -< A + « T - ,i. + 2(l-fl(l-oo.2rtD ' <26> 

Typically, « a < 1 ; f or example, = 0.001 and a = 0.2. The left- 
and right-hand sides of (26) as functions of 2w8T are sketched in Fig. 4. 
Apparently in the region of intersection, 2w8T <3C 13 and sin 2v8T 
« 2tST. Solving (26) with this approximation yields 

2tt(A + 8)Ttt^p8T. 

Thus 

.~ -ISA 



<* + 0' 



(27) 



and the necessary rate of rotation of the equalizer taps has been reduced 
by a factor of fi/(a + /3), which is about 1/200 for a typical case, 
a = 0.2, /3 = 0.001. The corresponding normalized excess mean squared 
error is 

If A = 1 Hz, a = 0.2, = 0.001, 7 1 = 1/2400 s. This amounts to 
about 10-*. 

V. STEADY-STATE SINUSOIDAL RESPONSE 

The phase jitter process {d„} that occurs in telephone channels is 
typically quasi-periodic. It is thus of interest to determine the steady- 
state solution of the coupled difference equations (13) and (14) when 
the driving term [6 n ] is sinusoidal. 

It is convenient at this point to rewrite eqs. (13) and (14) further 
in terms of eigenvalues and eigenvectors of the matrix a. Since Ct is 
Hermitian, its eigenvalues {ln}jL-w are positive real and its eigen- 
vectors (v.JfL-iv form an orthonormal set which is a basis in 2N + 1- 
dimensional space. Using these properties and expressing the vectors 
E„ and X as linear combinations of the {$<}, 

-, H 

i=-N 

X= £ Of*, 

i N 
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Fig. 4 — Illustration of the solution of 
2tt(A + S)T = 



-a/3 sin 2nST 



/3 2 + 2(1 -/3)(1 -cos27r5T) 

we can write (13) and (14) as 

S (B+ Di = (1 - 0\i)S ni e>6»+i-*»+ti + ^ ( e y(A n+l -A n+ i) _ i) 

A,- 



and 



A n+1 = a £ Im (O?,). 



- N -Zi^N (29) 
(30) 
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We now make the following change of variable in (29) and (30). 
Define 

ftfl^fc-W = Uni + jV ni . (31) 



Then we can write the real and imaginary parts of (29) as 

cos (0„+i - 

- N ^i^N (32) 



M (n+1)i = (1 - 0\i)Uni + ^~ [COS n ~ On) ~ COS (0„+l - 0„+i)l 



and 

\G\ 2 
v (n +i)i = (1 - Mvni + 1 y^ [sin (0„ - 6 n ) - sin (0 B+ i - n+1 )] 

-N^i^N, (33) 

and we can write (30) in the form 

N 

e„+i - dn = a £ [v ni cos (0„ - 0») - u ni sin (0« - 0»)J (34) 

i—N 

Equations (32), (33), and (34) are a set of nonlinear coupled differ- 
ence equations. In particular, eq. (34) is reminiscent of the equation 
governing a discrete-time, first-order, phase-locked loop. We shall solve 
linearized versions of (32), (33), and (34). Assuming the steady-state 
error angle (0„ — 0„) f or n ^> 1 is very small, we replace cos (6 n — 0„) 
by 1 and sin (0„ - 0„) by (0„ - 0„). Then (32) becomes 

U(n+l)i = (1 — {3\i)Uni, 

= (1 - p\i) n+1 u i, -N £ i ^ N, 

which approaches zero in the steady state (assuming /3 < 1/X, for all i). 
Thus in the steady state we are left with the linearized versions of (33) 
and (34) : 

V in +l)i = (1 - 0\i)Vni + ^-^- (0n+l -d n - n+ l + k) 

-N ^i^N (35) 
and 

d n +i - d n = a £ p ui . (36) 

i—N 

Equations (35) and (36) are linear and can be solved for a given 
sequence of channel phase shifts {0„}. We consider the case where the 
phase jitter is sinusoidal with frequency o> rad/s; i.e., 

d n = Re (Je*> nT ), 

where J is a complex constant. The solution for {v„.} is also sinusoidal : 

v^ = Re (Vie***): -N S i ^ N. (37) 
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Substitution of this trial solution in (35) and (36) yields a value of 
Vi after some algebraic manipulations. 

v = J{l-e^)\G^ 

X,(l - /3X,- - e* T ) (l - a £ IftlVQl - 0\* ~ e juT M^ 

(38) 

It follows from the sinusoidal variation of {v n i}iL-N that the error 
angle {0„ — 8 n \ and the equalizer tap coefficient vector C„ also vary 
sinusoidally with frequency o> in the steady state. 

The excess time-averaged mean squared error can be calculated from 
expression (31), (37), and (38). 

7 = < 7n > = <E;aE„) 

= E Ml«-*l"> 



\J\*\1 - e^ T \ l Si 
211 - e» T - aS 2 \ 2 ' 



(39) 



where 



and 



N \Gi\ 2 

Sl = ,.?* X,|l - p\ t /\l - e*r)\t 

N \GA 2 

The total mean squared error is, from (6b), 

<e„) = l - x*a-»x + 7 

n \n.\i |r|2ii_ pjoT\2a 
iJ^N \i + 2|1 -e* r - aS 2 | 2 k ' 

Typically, if the overall mean squared error is close to zero, 
Z J-^-'^l and 0Xi« |1 - «*"|. 

,--Ar A, 

Then the excess mean squared error in (40) is approximately 

\J\2\1 - e *-r|2 
2 1 1 — e» T - a\ 2 ' 

This expression corresponds to a previously derived, approximate, mean 
squared error due to sinusoidal jitter in the absence of noise [see eq. 
(39) of Ref . 3]. That equation, valid for a first-order, phase-locked loop, 
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was derived ignoring the coupling between eqs. (13) and (14) and 
assuming perfect equalization. Calculated curves of mean squared 
error versus a are found in Ref. 3. 

VI. CONCLUSIONS 

Previous studies have shown that the functions of joint passband 
equalization and data-directed carrier recovery in a qam receiver can 
be formulated as a gradient search algorithm. If the channel parameters 
entering into the expression for the gradient of the mean squared error 
are known, it is termed a deterministic gradient algorithm. In this 
paper we have analyzed the start-up behavior of the deterministic 
gradient algorithm and also the steady-state response to frequency 
offset and to sinusoidal phase jitter. The more practically motivated 
stochastic or estimated gradient algorithm, in which the channel 
parameters are initially unknown, has been studied experimentally and 
awaits further analytical study. 

It was shown that, under typical channel conditions, when the 
carrier phase offset is fixed, phase tracking does not greatly slow down 
the start-up behavior of the deterministic gradient algorithm, at least 
provided the equalizer adaptation coefficient /9 is much less than that 
of the phase estimator a. 

The phase estimator was first proposed as an adjunct to the pass- 
band equalizer, to mitigate the effects of too-rapid tap-coefficient 
rotation in the presence of channel frequency offset. It has been shown 
that frequency offset still causes tap rotation in the equalizer-plus- 
phase estimator system, but that the rate of rotation is tolerable, being 
on the order of 1/[1 + (a//3)] times the amount of frequency offset. 

The steady-state response of the linearized system to sinusoidal phase 
jitter was obtained. When linear distortion in the channel is not severe 
and the coefficient /3 is small, the system mean squared error due to 
tracking error approximates that of a first-order, phase-locked loop, 
as was assumed in an earlier paper. 

APPENDIX 

We wish to upper-bound the right-hand side of (16), given assump- 
tions (1), (2), and (3) of Section III. 

E; +1 ctE n+1 = e;(/ - /SCt)a(/ - /3a)E„ + x'a-^xie*--" - 1| 2 

+ 2 Re {E:(7 - /8a)X(l - «"*•«)}, (16) 

where A n+ i was given by (14). 

The first term on the right-hand side can be written 

E:aE„ - 20E* n a 2 E„ + /3 2 E;a 3 E„. 
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The matrix Ct is positive definite and Hermitian ; hence, 

-E:a 2 E n a -(a*E„)*a(a*E„) ^ -x min E;aE n , 

where X m i n is the minimum eigenvalue of a. Similarly, E^a 3 E„ 
^ XmaxE^aE„, where X max is the maximum eigenvalue. Thus we note 
for future reference that the first term in (16) is bounded as 

E* n (l - /3 a) a (7 - /3a)E„ ^ (l - 2/3X min + /3 2 x max )E:aE». (41) 

The second term in (16) is 

X*Q- 1 X|e^» +1 - II 2 ^ sin 2 



A„+i 



since X*a~ x X g 1. Upper-bounding sin 2 (A B+ i/2) by (A n+ i/2) 2 and 
substituting expression (14) for A n+ i, we have 

X*a- 1 X|e'*»+» - 1| 2 ^ a 2 [Im (E:X)] 2 . (42) 

The third term in (16) can be written as the sum of three terms. 

2 Re {E:(7 - 0a)X(l - erArt-i)} = 4 Re [Ej(/ - /3tt)X] sin 2 ^=±? 

- 2 Im (E*X) sin A n+ i + 2/3 Im (E*aX) sin A n+1 . (43) 

As in the inequality (42), the first term in (43) is upper-bounded by 

« 2 |E:(/-^a)X|[Im(E:X)] 2 . 

The matrix / — /3 a is Hermitian; its eigenvalues are {1 — /3X,-}, where 
the {X,-} are the eigenvalues of Cfc. Let X mB x and X m i n be the maximum 
and minimum eigenvalues, respectively. By assumption (3), 1 — /3X mtt x 
> and thus I — /3 a is positive definite. Therefore, 

|e;(/ - jsa)X| = |e;(/ - /3a)*a*a-*(7 - /3a)*x| 

^ [E;a*(7 - /3a)a*E B ]»[X*a-*(7 - /3a)a*x]*, (44) 

where we have used Schwartz's inequality. Using the positive definite- 
ness of I — /3 ft and a, we can further upper-bound the right-hand 
side of (44) by 

|E:(/ - /3Ct)X| g (1 - /3X min ) 2 (E:aE„)HX*a- 1 X)i 

^ (E:aE»)», (45) 

since the quantities 1 — /3X min and X*a _1 X are less than unity. Thus 
we have upper-bounded the first term in (43) by 

4 Re {E:(7 - PCL)X\ sin 2 ^ ^ a 2 (E;etE n )*[Im (E^)] 2 . (46) 
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After substituting for A n+ i using eq. (14), we can express the second 
term in (43) as 

-2 Im (E*X) sin A n+1 = -2a[Im (E^)] 2 sine [a Im (E*X)], (47) 

where 

sin 6 
sine 6 = — — • 
o 

The third term in (43) is 

2/3 Im (E^aX) sin A„ +i , 

which can be upper-bounded, using (14) and the inequality | sin A | 
^ I A | , by 

2a/3(|E:aX|)[|Im (E^X) | ] g e 2 /3 2 |E:aX| 2 

2 

-\ — j L~I m (E^X)] 2 for any arbitrary e, 
where we have used the simple inequality 

2apAB £ e 2 2 A 2 + - 2 B\ 
But 

|E:ax| 2 = |(E;a*)(a-*X)| 2 
^ E* n a 3 E n , 

by Schwartz's inequality and the fact that X*Ct -1 X ^ 1. 
Thus the third term in (43) is upper-bounded by 

e 2 /5 2 (E;a 3 E„) + 2j [Im (ESX)] 2 . (48) 

Finally, substituting (42), (46), (47), and (48) into the right-hand 
side of (16), we have 

7„+i = E; +1 aE n+1 ^ E:aE n - 2/?E:a 2 E„ + /3 2 (i + e 2 )E:a 3 E„ 

+ a(R„[Im (E:X)] 2 , (49) 
where e is arbitrary and 

Gt» = a \l + (E^GEn)* + j 2 1 - 2 sine [a Im (E^X)]. (50) 

We make the following choice of e : e = eo, where eo is defined by (with 
70 = EottEo) 

a(l + ^ + - 2 \ = 2 sine («V^). (51) 
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Figure 2 is a sketch of the left- and right-hand sides of eq. (51) as func- 
tions of a for various values of eo- Equation (51) has a unique solution 
with ^ eg < °° as long as ^ a < a , where a is denned by 

ao(l + V70) = 2 sine (oioVyo). 

Note also that, by assumption (3), the coefficient 1 — 2/3X m j n 
+ 2 Xmax of E*aE„ in the bound (41) is less than 1 and hence (49) 
can be weakened to 

E; +1 aE„ +1 ^ KaE n + <R„[Im (EiJX)] 2 . (52) 

Lemma: (R„ is negative, and hence the sequence {y n = E* n dE n ) is mono- 
tone decreasing. 

Proof: We first observe that the sine function in (50) defining (R n is 
even, positive, and monotone decreasing provided its argument's 
absolute value is less than w. But its argument is 

aIm(E»X) ^ a|Eja*a-*X|. 

This can be bounded, using Schwartz's inequality, by 

a (E* n aE n x*a- l X)* ^ a (E*aE„)i 

and so 

- sine [a Im (EjX)] ^ -sine |>(E*aE»)] for a(E:aE n )* < t. 

(53) 

In particular, 

a Im (E$X) ^ aVr^ < tt 

by assumption (1), and hence we can upper-bound (Ro by 

(Ro £ a ( 1 + Vyo + p J - 2 sine (ajyo). (54) 

According to our choice of e = eo, defined by (51), the right-hand 
side of (54) is zero, and so (R ^ 0. It follows from (52) that V71 ^ V70, 
which is less than r by hypothesis. Thus (Ri is bounded, using (53) and 
e = 60, by (Ri ^ (Ri, where (R„ is defined by 

(R„ = <* Tl + 7 | + 1 1 - 2 sine (orS), (55) 

and 01 = by the definition of el. Now since V71 ^ V70 ^ t, 
— 2 sine (ay\) < —2 sine {ay$) 

and so 

Bi ^ & ^ ^0 = 0. 
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Similarly, from (52), 72 < 71 and by induction 

7n ^ 7n-i ^ • • • ^ 7o 
and all (R„ ^ 0. 



Q.E.D. 



Finally, since (R„ is negative, we obtain the following recursive upper 
bound from (49) : 

7 n+i = E; +1 aE n+1 ^ E* n aE n - 2pE* n am n + /3 2 (1 + e§)E;a 3 E„. (56) 
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